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Abstract. We report the results from a spatially resolved spectroscopy study with XMM-Newton on the relaxed cluster of 
galaxies Abell 478. From the EPIC data we extract a temperature profile and radial abundance profiles for Ne, Mg, Si, S, Ca, Fe 
and Ni. The abundance profiles follow the same trends as observed in other clusters. The spectra of the core of the cluster can 
be best fitted with a multi-temperature model. We argue that this multi-temperature behavior is mostly due to projection effects, 
because of the strong temperature gradient in the core. Contributions from other effects, for example, intrinsic temperature 
stratification cannot be fully excluded. For the first time we measure an underabundance of oxygen in the Galactic absorption 
component towa rd a cluster. The measured oxygen abundance in this absorber is about 0.5 times the solar oxygen abundance 
as determined bv lAnders & Grevessd I\98& . 

Key words. Galaxies: clusters: general - Galaxies: clusters: individual: Abell 478 - cooling flows - intergalactic medium - 
X-rays: galaxies: clusters - ISM: abundances 



I. Introduction 

The evolution and hydrodynamical structure of the hot diffuse 
X-ray emitting gas in clusters of galaxies is still not well un- 
derstood. Several clusters show excess X-ray emission in their 
cores. The s o-called cooling flow models (for a review see 
lFabianlfl994l) propose radiative cooling of gas in the core to 
explain the observed surface brightness and temperature pro- 
files. The pressure decrease in the core due to the cooling gas 
causes a net inflow toward the center of the cluster, hence the 
name cooling flow. XMM-Newton observations show that the 
amount of cool gas in the core of cooling-flow clusters is much 
lower than predicted jPetolcmjM^^ 

iTamura et all2001bllPeterson et all2003tlKaastra et all2004 . 
Therefore, the standard cooling-flow model needs to be ad- 
justed. Several mec hanisms to explain t he observations have 
been proposed (see IPeterson et alJl200li for a list of possible 
explanations). 

The Abell 478 cluster of galaxies is a good exam- 
ple of a highly relaxed cluster. Earli er X-ray observations 
with EXOSAT fedge & Stewart! 1 199 lb. Ginga and Einstein 

II. Tohnstoneetalll992h.ROSATl|Al 
11994 , ASCA jMarkevitch et al 



en el 
1998 



lJl993|lwhite et alJ 
IWhitd 12000) and 

Chandra ( Su n et alJl2003l) show a significant excess of X-ray 
emission in the core, suggesting the presence of a massive cool- 
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ing flow. The observations performed by missions with suffi- 
cient spatial resolution (e.g. ROSAT, ASCA and Chandra) are 
consistent with a radial temperature decrease toward the core. 
In a recent Chandra observation an X-ray cavity was discov- 
ered in the core which se ems to be associated with a lobed radio 
source JSun et al.l2003l) . 

In this paper we study the properties of Abell 478 using 
high-resolution and spatially-resolved spectra o btained with 
the E uropean Photon Imaging Camera (EPIC, iTurner et alJ 
2001) and the Reflection Grating Spectrometer (RGS, 
Iden Herder etalJ 1200 ll) aboard XMM-Newton (Jansenetal. 
2001J). We focus on the physics and metal abu ndances in the 
core o f Abell 478. In an associated paper Poin tecouteau et all 
(2004) describe the large scale temperature, gas and dark mat- 
ter distribution in this cluster using a pointed and an offset 
XMM-Newton EPIC observation. 

Throughout this paper we use Ho = 50 km s _1 Mpc~' and qo 
= 0.5. Using this cosmology 1' is 157 kpc at the cluster redshift 
of 0.0881. 



2. Observations and data analysis 

The observation of Abell 478 was performed as part of the 
Guaranteed Time program on February 15 2002 and had a to- 
tal duration of 126 ks. Both EPIC MOS cameras were operated 
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Table 1. Maximum count rates and resulting useful exposure 
times after screening for solar flares. 



Detector Max. count rate Useful exposure time 

(counts s -1 ) (ks) 
PN 053 51 

MOS 51 
RGS 0.36 119 



Table 2. Systemat ic errors in cluded in the EPIC spectral fits. 
(Adapted from Kaa stra et al]2 004) 



Energy band Error source Error background 
0.4-0.5 10% 35% 
0.5-0.7 5% 25% 
0.7-2.0 5% 15% 
2.0-10.0 5% 10% 



in Full Frame mode and the EPIC pn camera in Extended Full 
Frame mode. For all EPIC cameras the thin filter was used. 

For the analysis we use the event lists produced by the 
XMM-Newton SOC and the 5.4. 1 version of the XMM Science 
Analysis System (SAS). To correct for enhanced and variable 
background due to soft proton flares we use the events with 
energies larger than 10 keV from pn and CCD 9 of RGS. The 
data are filtered using upper and lower count-rate thresholds. 
The threshold s are calcu l ated in a way analogous to the method 
used in IPratt & Arnaudl J2002h . We first make a histogram of 
the light curve above 10 keV with a binsize of 100 s and pat- 
tern equal to 0, which we fit it with a Poissonian function. The 
count rate thresholds are subsequently fixed to A^± 3 V^V, where 
N is equal to the mean number of counts within a bin. Because 
the pn is most sensitive to solar flares, we merge the Good 
Time Intervals (GTI) derived from the pn data also with the pre- 
liminary GTI for MOS. Therefore pn and MOS have roughly 
the same exposure time. The maximum count rates and the re- 
sulting useful exposure times are shown in Tabled This sim- 
ple variable background procedure aims at obtaining maximum 
statistics in the core region, where the background is relatively 
less important, but excludes contributions from the large solar 
flares. Because the cluster signal in the 0^4' region dominates 
over the background, it is not necessary to perform a double 
background subtra ction nor a multi- band flare search like de- 
scribed in lPointecouteau et alJJ2004l) . The countrate thresholds 
we chose are well within the advised thres hold of < 1.0 c ts s _1 
(PN) provided by the XMM-Newton SOC jKirschl2003l) . 

2.1. EPIC 

Because the cluster is slightly elongated along the north-east 
to south-west axis, we extract the source and background spec- 
tra from elliptical annuli allowing patterns <12 for MOS and 
<4 for pn. We use a ratio of 1.2 between the major and mi- 
nor axis of the ellipse, which we obtain from an empirical 2D 
/?-model fit to the EPIC image. Although the data show some 
minor deviations fro m this 2D model, the value is in the range 
of 1.2-1.4 (RO S AT. IWhiteetal.ll 1994 an d compatible with 
1 .22 derived by Poin tecouteau et al.l (120041) who uses the same 
method and data. Further on in this paper we characterize the 
extraction region by its semi major axis. We extract the EPIC 
backgrou nd spec t ra from the blank sky event files provided by 
iRead & Ponmanl J2003). Relatively bright point sources in the 
EPIC field of view have been excluded both from the Abell 478 
and background datasets. 



To gain statistics in each spatial bin we extract the spec- 
tra from annuli with a width of more than 30". This way, we 
are less sensitive to the energy dependent shape of the PSF and 
therefore we neglect this effect in the rest of the analysis. The 
downside of choosing large binsizes is that mixing of temper- 
atures within the bin also increases with binsize if there are 
strong temperature gradients present. This effect may become 
important during multi-temperature fitting. 

Although we neglect the influence of the energy dependent 
part of the PSF, a substantial broadening effect in the core re- 
mains. To correct the normalization of each annulus for PSF 
effects we calculate correction factors by fitting the surface 
brightness profile from the Advanced Camera Imaging System 
(ACIS) on board Chandra dSun et all2003l) with a double King 
profile. We normalize both the ACIS fit and the EPIC surface 
brightness profile and calculate the ratios between the two for 
each bin. These factors range from 0.73 in the core up to 1.03 
at 4'. Then, we multiply the effective area of each annulus with 
this factor and obtain a PSF corrected response. 

For the spectra l analysis we use the SPEX package 
jKaastra et alJl2003h and fit the spectra over the 0.4-10 keV 
range. We include systematic errors due to uncertainties in the 
calibration and the background in the spectral fit. The applied 
syste matic errors ar e shown in Table|2] they have been adapted 
from iKaastra et alJ ll2004l) . These systematic errors enclose 
possible variations of the Cosmic X-ray Background (CXB), 
background normalization and errors in the calibration. 

After a preliminary analysis, the pn appears to have a gain 
problem, resulting in a significantly lower value for the red- 
shift of the source than measured by the MOS instrument. From 
MOS we obtain a redshift of 0.0889 + 0.0003, which is consis- 
tent w ith the optical value of 0.0881 ± 0.0009 dZabludoff et alJ 
1990). But, our fits of the pn spectra result in a redshift of 
0.0775 ± 0.0002. To correct for this, we measure the centroid 
energies of the aluminum, nickel, copper and zinc background 
lines. A linear fit (E nev , = a£oid + b) through the measured and 
expected energies for these lines yields a — 1 .0066 ± 0.0007 
and b = (-1.04 + 0.24) x 10~ 2 . After correcting the energy val- 
ues in the original eventfile according to this relation, we fit the 
resulting spectra again. The corrected spectra show large dis- 
crepancies at low energies. Therefore, we conclude that the ob- 
served gain pro blem is probably not linear , in line with the ar- 
gumentation of lPointecouteau et alJ J2004). Unfortunately, the 
number and distribution of background lines is insufficient to 
derive the exact shape of the gain correction. For the rest of the 
analysis we therefore use the original event file, fit MOS and 
pn separately and let the redshift free in the fit. 
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Fig. 1. Single-temperature fit results. EPIC MOS (stars), pn (diamonds) and RGS (plus-sign) are shown. The dashed line shows 
an empirical powerlaw fit to the Fe abundance distribution for comparison. 



2.2. RGS 

We extract the RGS spectra with SAS version 5.4.1 following 
the same method as described in Tam ura et alJ J2001al) . We se- 
lect the events from a rectangular area on the CCD strip with 
a full width of 1' in the cross-dispersion direction and centered 
on the core of the source. We extract t he background spectra 
from several blank field observations ( Tamur a & den Herded 
120031) . Because the observation was performed before the 
cooling of the RGS instruments, the spectra contain a rela- 
tively high number of bad columns due to warm pixels. These 
bad columns appear as absorption-like features in the spectra. 
Unfortunately, also CCD gaps are present in the wavelength 
ranges where line emission is expected. Together with the rela- 
tively high thermal Bremsstrahlung component due to the high 
temperature of the cluster, the lines are difficult to resolve. 
Therefore, we also include the second order spectrum in the 
fit to increase the statistics. All RGS spectra are fitted over the 
8-35 A range. 

Because the RGS gratings operate without a slit, the result- 
ing spectrum of an extended source is the sum of all spectra 
in the (in our case) l'x ~ 12' field of view, convolved with 
the PSF (see for a complete discussion about grating responses 



lDavisll200ll) . Extended line-emission appears to be broadened 
depending on the spatial extent of the source along the disper- 
sion direction. In order to describe the data properly, the spec- 
tral fits need to account for this effect. In practice, this is ac- 
complished by convolving the spectral models with the surface 
brightnes s profile of the source along the dispersion direction 
( Tamu raet alJ 12004 . For that purpose we extract the cluster 
intensity profile from MOS 1 along the dispersion direction of 
RGS, which we convolve with the RGS response during spec- 
tral fitting. Because the radial profile of an ion can be different 
from the mean profile, this method is not ideal. We let the scale 
of the width and the position of the profile free in the fit to 
match the profiles of the main emission lines. 



2.3. Spectral models 



We fit the spectra both with a single temperature collisionally 
ionized plasma model (MEKAL) and a differential emission 
measure (DEM) model called wdem. In this particular model 
the emission measure (Y) is distributed as a function of tem- 
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Table 3. Averaged EPIC results for Abell 478 for the single temperature (1) and wdem (2) fits. Errors are given at the Icr 
confidence level and the <-sign denotes a 2<x upper limit. We leave the abundance of oxygen in the absorption component free to 
fit the oxygen edge near 0.5 k eV. This parameter is called <9 a b s - The abundances are given relative to solar abundances measured 
bv lAnders & Grevessd(ll989l) . 





Model 






EPIC 










0-0.5' 


0.5-1.0' 


1.0-2.0' 


2.0-3.0' 


3.0^1.0' 


X 2 /d.o.f. (MOS) 


1 


318/327 


286 / 327 


290 / 327 


288 / 327 


281 /303 




2 


247 / 326 


285 / 326 


292 / 326 


278 / 326 


281 /302 


v 2 /H o f ( nnl 


1 


216 / 179 


204 / 179 


228 / 179 


211 / 179 


166 / 158 




2 


177 / 178 


193 / 178 


228 / 178 


209 / 178 


161 / 157 


N H (10 21 cm -2 ) 


1 


3.54 ± 0.06 


3.40 ± 0.05 


3.32 ± 0.03 


3.34 ± 0.06 


3.49 ±0.11 




2 


3.60 ± 0.04 


3.44 ± 0.04 


3.32 ± 0.04 


3.34 ± 0.05 


3.52 ± 0.05 


O u 


1 


0.40 ± 0.02 


44 + 03 


0.41 ± 0.02 


39 + o 04 


32 + 05 




2 


0.50 ± 0.02 


0.45 ± 0.02 


0.43 ± 0.03 


0.38 ± 0.03 


0.34 ± 0.06 


kT (keV) 


1 


4.70 ± 0.03 


5.88 ± 0.04 


6 30 + 04 


6 76 + 07 


6.42 ± 0.07 


kT m „ (keV) 

' v * max \ lvw * / 


2 


8.31 ± 0.08 


7.7 ± 0.3 


7.0 ± 0.4 


8.8 ± 0.4 


9.5 ± 0.5 


1 /a 


2 


i.iz _ 03 


n 3Q + n 1 5 


1? + 03 


40 + 1 1 


0.8 + 0.4 





1 


< 0.06 


< 0.14 


< 0.09 


< 0.09 


< 0.11 




2 


< 0.07 


< 0.06 


< 0.09 


< 0.06 


< 0.05 


Ne 


1 


1.32 ± 0.11 


0.90 ± 0.12 


0.63 ±0.11 


< 0.5 


1.04 ± 0.17 




2 


< 0.18 


0.52 ±0.11 


60 + 1 1 


0.4 ± 0.2 


43 + 016 


Ms 


1 


0.49 ± 0.19 


0.63 ± 0.15 


0.51 ± 0.15 


< 0.5 


0.53 ± 0.19 




2 


< 0.08 


0.46 ± 0.14 


0.49 ± 0.15 


< 0.4 


0.31 ± 0.17 


Si 


1 


0.23 ± 0.07 


0.40 ± 0.08 


0.20 ± 0.07 


< 0.5 


0.36 ± 0.16 




2 


0.22 ± 0.06 


0.35 ± 0.07 


0.20 ± 0.07 


< 0.4 


0.32 ± 0.13 


s 


1 


0.38 ± 0.08 


0.25 ± 0.09 


0.22 ± 0.09 


0.24 ± 0.12 


<0.08 




2 


0.41 ± 0.06 


0.25 ± 0.09 


0.22 ± 0.09 


0.23 ±0.11 


< 0.09 


Ar 


1 


<0.4 


< 0.2 


<0.2 


<0.15 


< 0.6 




2 


0.35 ±0.15 


<0.2 


<0.2 


<0.16 


< 0.6 


Ca 


1 


0.73 ±0.19 


0.8 ± 0.2 


0.7 ± 0.2 


< 0.8 


1.1 ± 0.4 




2 


1.1 ± 0.2 


1.0 ± 0.2 


0.7 ± 0.2 


< 0.9 


1.2 ±0.4 


Fe 


1 


0.410 ± 0.010 


0.351 ± 0.009 


0.320 ± 0.008 


0.285 ±0.011 


0.268 ± 0.014 




2 


0.425 ± 0.010 


0.367 ± 0.009 


0.322 ± 0.008 


0.294 ±0.011 


0.292 ± 0.014 


Ni 


1 


1.2 ±0.2 


0.78 ± 0.19 


0.55 ±0.17 


< 0.8 


0.9 ± 0.3 




2 


0.79 ±0.19 


0.7 ± 0.2 


0.54 ± 0.17 


< 0.8 


0.7 ± 0.3 



3. Results 

3.1. EPIC 

The results of the EPIC single-temperature fits are presented in 
Fig-E an d Table [5] The upper left plot in Fig^shows that the 
absorption is not a constant across the cluster. The temperature 
profile shows a steep decrease toward the core with a hint for 
a bend outside the 2.0' radius. The abundances of neon, sul- 
fur and iron are consistent with an increase toward the center. 
However, neon and sulfur are also consistent with a flat distri- 
bution, like the other abundances. Because of low statistics we 
had to discard the 8-10 keV band in the 3-4' annulus. 

Fig- 12 shows the residuals of the combined MOS spectra 
of the 0-0.5' region, fitted with a single temperature compo- 
nent for which the abundances have been set to zero afterward. 
In this way we show the lines that are detected in the EPIC 
spectrum. Except for neon, the abundance of all elements can 
be fitted independently. Since the main neon line is close to 
the Fe-L complex blend at ~1 keV, it is difficult to constrain 
the neon abundance with EPIC. Because the Fe-K line is the 
strongest line, this line mainly determines the iron abundance. 



perature (T) as shown in Eq. (fl) adapted from iKaastra et all 

fcnoA : 

dY = (cT a r<r max 

dT \ o t > r max { ' 

This model is an empirical parametrization of the DEM dis- 
tribution found in the core of many clusters. In this form the 
limit a — > oo yields the isothermal model. For convenience we 
will use 1/ff in this paper, because then the isothermal model 
is obtained when I /a = 0. The classical cooling-flow model 
corresponds to a dY/dT of 1/A(T), where A(T) is the cooling 
function. 

To investigate the role of projection effects, we fit the spec- 
tra of the core also with an extra temperature component of 6.5 
keV. This way we emulate the spectral contribution of the outer 
parts of the cluster. 

We notice from a preliminary analysis that the oxygen edge 
near 0.5 keV is not well fitted. Therefore, we use an absorp- 
tion model component with variable element abundances and 
leave the oxygen abundance free. We call this particular oxy- 
gen abundance O a bs further on in the paper. 



J. de Plaa et al.: X-ray spectroscopy on Abell 478 with XMM-Newton 



5 



Table 4. Averaged results from the single-temperature (3) and 
wdem (4) fits with an extra 6.5 keV temperature component. 





Model 


EPIC 






0-0.5' 


0.5-1.0' 


^ 2 /d.o.f. (MOS) 


3 


244 / 326 


284 / 326 




4 


247 / 325 


282 / 325 


^ 2 /d.o.f. (pn) 


3 


178 /178 


192 / 178 




4 


176 / 177 


200 / 177 


Nh (10 21 cm- 2 ) 


3 


3.60 ± 0.03 


3.45 ± 0.06 




4 


3.58 ± 0.03 


3.43 ± 0.06 


o abs 


3 


0.453 ± 0.016 


0.41 ± 0.04 




4 


0.488 ± 0.017 


0.48 ± 0.03 


kT (keV) 


3 


2.46 ± 0.04 


4.14 ±0.08 


kT mia (keV) 


4 


3.04 ± 0.05 


4.95 ± 0.17 


1/ar 


4 


0.38 ± 0.07 


1.0 ± 0.3 


O 


3 


< 0.08 


< 0.12 




4 


< 0.08 


< 0.05 


Ne 


3 


83 + 1 6 


65 + 14 




4 


0.26 ± 0.12 


0.20 ± 0.10 


Mg 


3 


77 + 07 


50 + 1 5 




4 


< 0.09 


0.43 ± 0.13 


Si 


3 


23 + 05 


36 + 08 




4 


0.22 ± 0.05 


0.33 ± 0.06 


S 


3 


38 + 06 


27 + 09 

VJ • 1 _!_ \J . \J J 




4 


0.39 ± 0.06 


0.26 ± 0.08 


Ar 


3 


0.31 ±0.14 


< 0.18 




4 


0.34 ±0.15 


<0.2 


Ca 


3 


1.1 ± 0.2 


1.0 ±0.2 




4 


1.09 ±0.19 


1.0 ± 0.2 


Fe 


3 


0.436 ± 0.010 


0.362 ± 0.009 




4 


0.429 ± 0.010 


0.363 ± 0.009 


Ni 


3 


0.83 ±0.18 


0.7 ± 0.2 




4 


0.77 ±0.19 


0.69 ±0.19 



The strength of the Fe-L complex is mostly dependent on the 
temperature structure. 

In the core of the cluster the x 2 value from the single- 
temperature fit is higher than can be expected from the high 
signal-to-noise in that area. We therefore also fit the core with 
the multi-temperature wdem component defined in Sect. 12.31 
These results are also listed in Table [3] The wdem abundances 
are all consistent with the single temperature abundances, apart 
from neon, which is blended by Fe-L. From the x 2 values 
the multi-temperature model appears to be preferable over the 
single-temperature model in the 0-0.5' region. 

Then, we add a single temperature component fixed at 6.5 
keV to emulate the projection effect from the hot gas in front 
of the core. The value of 6.5 keV is roughly the mean temper- 
ature outside the 1' radius. In Table |4] we show the results of 
this fit. The abundances are again consistent apart from neon. 
We notice that the fitted temperatures for the core are now sig- 
nificantly cooler than without the hot component. The value for 
1/a in the wdem component is 0.38 + 0.07, lower than the value 
we previously found, 1.12+Q3 2 . 



Table 5. Results from the RGS spectral fits of the core re- 
gion (0-0.5') of Abell 478. The RGS spectra were fitted using 
the wdem model component and a single-temperature model. 
Errors are given at the l<x confidence level and the <-sign de- 
notes a 2cr upper limit. We left the abundance of oxygen in the 
absorption component free to fit the oxygen edge near 0.5 keV. 
This parameter is called O a b s . 





RGS 


RGS 


Unit 




wdem 


single-temperature 




X 2 /d.o.f 


867 / 753 


887 / 754 




Nh 


3.15 ±0.08 


3.30 ± 0.07 


10 21 cm" 2 


o abs 


0.67 ± 0.04 


0.63 ± 0.04 


z© 


kT 




3.5 ± 0.4 


keV 


kT 

max 


11 r 12 

11.1 23 




keV 


1/a 


1.3 ±0.3 






O 


0.15 ± 0.07 


0.09 ± 0.05 


Z 


Ne 


0.46 ±0.13 


0.39 ±0.10 


z 


Mg 


0.17 ± 0.12 


<0.2 


Zq 


Fe 


6r 010 


0.23 ± 0.06 


z 



3.2. RGS 

We obtain a reasonable fit to the RGS spectra of the core (0- 
0.5') with the single-temperature model. The x 2 I d.o.f. is 887 
/ 754. Unfortunately, the O vra line complex is partly within a 
CCD bad-column on RGS 1 and falls in the dead area of RGS2, 
which makes the derived oxygen abundance more uncertain. 
The Fe-L complex and Ne line near 12 A are detected and re- 
solved. 

Because the EPIC spectra from the core region are best fit- 
ted with a multi-temperature model, we also fit the RGS spec- 
tra with the wdem component. This results in a slightly better 
X 2 I d.o.f. of 867 / 753. The results of this fit and the single- 
temperature model are shown in Table [5] and Fig. [3] Because 
the temperature determination with RGS depends mostly on the 
line emission, the weakness of the lines in this spectrum causes 
the temperature to be less well constrained. Furthermore, the 
energy range of the RGS is much smaller than EPIC and 
concentrated on the soft energy part of the spectrum, mak- 
ing it more sensitive to the cool component of the spectrum. 
Therefore, the fit is less sensitive to the multi-temperature 
distribution. Because of the increased freedom of the fit due 
to the use of multi-temperature components, the wdem fit is 
marginally better than the single-temperature fit. We have also 
attempted to fit the RGS spectrum with an extra hot component 
of 6.5 keV which we added to both a single-temperature and 
wdem model. Because the temperatures in this model could not 
be constrained, we have not pursued this model any further. 

4. Discussion 

We analyze the high-resolution spectra of the Abell 478 cluster 
of galaxies and we derive, for the first time, radial abundance 
profiles for several elements in this hot cluster. The exposure 
time of 126 ks allows us to resolve the main emission lines. 
By fitting these lines and the continuum emission we are able 
to put constraints on the temperature structure in the core and 
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Fig. 2. Residuals of the averaged MOS spectrum of the core of Abell 478 with the abundances put to zero. 



the abundance distribution of Ne, Mg, Si, S, Ca, Fe and Ni. 
Furthermore, by fitting the oxygen edge near 0.5 keV we mea- 
sure an underabundance of oxygen in the absorption compo- 
nent. 

The temperature profile obtained from the single- 
temperature fits (see Fig.^ an d Table|3 of the core (< 2') are 
consistent with th e temperature profile measured with Cha ndra 
llSun et alJl200l and ROSAT fWhite etaljfl994l lAllen et alJ 
Il993l) . It confirms that Abell 478 has one of the steepest tem- 
perature gradients observed so far. 

In the outer regions, outside 2' from the core, our pn 
temperature profile shows a hint of a bend which is con- 
sistent with the more exte nded temperature determinati on by 
IPointecouteau et all d2004l) . but is inconsistent with ISun et alJ 
J2003h . The ACIS aboard Chandra, which has a lower sensi- 
tivity at high energies and a smaller field of view than pn, is 
therefore more vulnerable to systematic background effects in 
the temperature. The bend in the pn temperature profile is com- 
parable with profiles ob served in other clus ters, for example the 
core of Sersic 159-03 JKaastra et all200ll) . and has been inter- 
preted as the transition from the cooling core to the inter-cluster 
medium. 

In the core region, within 0.5', a multi-temperature (wdem) 
model fits better than the single-temperature model. The value 
of 1/a — 1.3 ± 0.3 from RGS is high compared to the cluster 



sample s studied by e.g. lPeterson et al.(j2003b and lKaastra et alJ 
(2004). These studies find values for 1/a which are roughly of 
the order of 0.5. If we add an extra temperature component 
of 6.5 keV to account for the emission from the outer regions 
in front of the core, then the difference between the single- 
temperature and wdem component vanishes (see Table0}. This 
suggests that the projection effect of the hot cluster material 
along the line of sight strongly affects the 1/a value in the fits. 

The multi-temperature behavior is not necessarily ex- 
plained by contamination from the oute r parts of t he clu ster. 
Although the X-ray cavities reported bv lSun et alJ J2003h are 
not spatially resolved with XMM-Newton, the value of 1/a and 
the cool area associated with the X-ray cavity may be related. 
Together with the high resolution study of the temperature pro- 
file (IPointecouteau et a"fll2004l) . which reveals a steep temper- 
ature gradient in the core region, these effects could also par- 
tially account for the multi-temperature behavior observed in 
the central bin. Unfortunately, the spatial resolution of XMM- 
Newton does not allow us to draw conclusions from these data 
on the complex structure in the core. 

From the single-temperature fits we measure a central in- 
crease of the iron abundance from ~0.3 at about 4' to ~0.4 in 
the center which is consi stent with the mean abundance profile 
of Abell 478 derived bv lWhitd §000) using ASCA. The cen- 
tral abundance is also consistent with earlier EXOSAT mea- 
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Fig. 3. RGS averaged spectrum (0-0.5') of Abell 478. Only the first order spectrum is shown for clarity. Data points from bad 
columns were removed from the RGS 1 and RGS2 spectra before they were averaged. 



surements by Edg e & Stewart! Jl99ll) . The significant central 
increase of Fe is similar to the profiles of other hot (> 6 keV) 
clusters ( Tam ura et all200 4). Also neon and sulfur seem to fol- 
low the same trend in our results. However, due to the large 
error bars these profiles could still be consistent with a flat dis- 
tribution. The cent ral increase o f sulfur was also found in the 
cluster sample of Tam ura et alJ |.2004). Neon is not resolved 
in EPIC, because of the lower spectral resolution and blending 
with lines from the iron L complex. Therefore, the systematic 
uncertainty on the abundance of neon is quite large. Moreover, 
the single-temperature model tries to fit the Fe-L complex by 
enhancing the Ne abundance, while the reason for the Fe-L to 
be enhanced is the presence of colder material. Because RGS 
resolves the neon line, its value for the neon abundance is much 
more robust. Neon and oxygen are thought to have the same 
origin. Therefore, we expect their values and spatial behavior 
to be comparable, which excludes the very high numbers for 
neon from the single-temperature EPIC fits. The nickel abun- 
dance we find is likely to be overestimated du e to an error in 
the nickel line energies (see also lGastaldello & Molendil20 04). 

The oxygen abundance could not be constrained in our 
EPIC fits, but despite the occurrence of bad columns in the 
line resolved by the RGS we could get a value for the oxy- 



gen abundance consistent with the EPIC upper limits. We can 
compare our result with the predicted O/Fe ratios from the- 
oretical models for supernovae la and II. The O/Fe ratio ob- 
tained from the RGS wdem fit is 0.25 + 0.13, which is lower 
than the average va lue of 1.2 ± 1.6 for hot clusters reported by 
iTamura et alJ (2004). but not significantly different because of 
the large spread. The theoretical models predict an O/Fe ratio 
of < 0.05 and 1 .5 — 4 for supernova types la and II, respectively. 
Our value of 0.25 + 0. 13 is in between the two predictions, like 
in the other hot clusters. 

In Table [3] and Table [5] the values for A^h and O a b s appear 
to be different between EPIC and RGS. This is partly due to 
the fact that the EPIC/RGS cross-calibration near the oxygen 
edge is still not optimal. The EPIC results are also slightly af- 
fected by correlations between A^h, O a b s and the oxygen abun- 
dance, because of blending due to the lower CCD resolution. 
The O a bs anti-correlates with absorption, but because fixing 
the edge results in large residuals at lower energies, we de- 
cided to let the oxygen abundance in the absorption component 
free. Fig.0]shows the error ellipses for the oxygen abundance 
against Nh- From the plot it is apparent that Nh and O a b s are 
anti-correlated. If we fix the O a b s abundanc e to 1 .0, the Nh val- 
ues we find are consistent with those from Poin tecouteau et all 
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fc0u4 and lSun et aD d2003l) . T his plot also shows that O abs is 
significantly different from solar. Weisskopf et alJfc004l) report 
from their Chandra obs ervation of the Crab nebul a that, us- 
ing the abundances from lAnders & Grevessd {l989), they find 
an underabundance of o xygen in the Galac tic absorption com- 
ponent of 0.41 ± 0.07. iTakei et alJ fc003l) find values in the 
range from 0.63 ± 0.12 to 0.74 ± 0.14 using a Chandra ob- 
servation of Cyg X-2. Within the calibration uncertainties our 
results are consistent with these numbers. It is likely that the 
observed underabundan ce is an art ifact of the solar abundances 
we use. lAllende Prieto et alJ J200lh reported a new solar photo- 
spheric abundance of log e(O) = 8. 69 ±0.05 dex. This is 0.58 + 
0.08 times smaller than the value o f Anders & Grevesse ( 1989) 
and in between our values for O a b s . Pointec outeau et al. ( 2004) 
show that the origin of the absorption in the direction of Abell 
478 is most likely Galactic and correlated with infra-red data. 
A lthough O a hs is consistent with the solar oxygen abundance 
of Ellende Pri eto et alJ fcOOlb we cannot fully exclude the ex- 
istence of a real underabundance of oxygen in the absorbing 
material along the line of sight. 

Finally, we find that the average redshift measured in 
the pn detector is 0.0775 ± 0.0002 compared to 0881 + 
0.0009 measured in the optical by Zabludo ffet alJ ([1990). 
This is probably a gain-related problem which was not cor- 
rected for during the XMM-Newton SOC data processing. 
Problems with the redshift determinations are not unique. 
Discrepancies between X-ray measur ed redshifts and op tical 
redshifts have been reported ea rlier bv lZhang et alJ J2004I) and 
iTakahashi & Y amashita J2003I) . 



5. Conclusions 

We analyze the high-resolution XMM-Newton spectra of the 
Abell 478 clusters of galaxies and conclude that: 

- We measure an underabundance of oxygen in the Galactic 
absorption component. This is measured for the first time 
using an observation of a cluster of galaxies. 

- We derive radial abundance profiles for Ne, Mg, Si, S, Ca, 
Fe and Ni , which con firm the trends observed in other clus- 
ters JTamura et all2004l) . 

- The core of the cluster shows multi-temperature behavior, 
which is mostly explained by projection effects partly be- 
cause of the steep temperature gradient. 
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